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Abstract 

The calculation of the transport coefficients of the strong quark-gluon plasma 
' r-| '- (sQGP), i.e. sheer viscosity t], heat conductivity k and bulk viscosity (, in 

Q-1 the first Chapman-Enskog approximation is presented. Their formulation 

in terms of two-fold integrals depending on the particle interaction, known 
as relativistic omega integrals, is derived and their evaluation in the ultra- 
relativistic limit is worked out assuming a cross section independent of the 
relative and total momentum of two colliding particles. We find a suppression 
of the bulk viscosity and a pronounced temperature dependence for the sheer 
viscosity and heat conductivity. However, at high temperatures, they scale 
with the third and second power of the temperature respectively as expected. 
Furthermore, we find that all results in this ultrarelativistic expansion are 
dominated by the leading order contribution. 

Key words: Quark-gluon plasma. Viscosity, diffusion, and thermal 
conductivity. Kinetic and transport theory of gases 
PACS: 12.38.Mh, 51.20.+d, 51.10.+y 



1. Introduction 



Transport properties of the strongly interacting quark-gluon plasma (sQGP), 
in particular the bulk and shear viscosities, which describe the hydrodynamic 
response of the system to energy and momentum density fluctuations, have 
attracted remarkable attention in the last years. The experimental findings of 
the heavy-ion collisions at the Relativistic Ion Collider (RHIC) have pointed 
to a small specific sheer viscosity, i.e. the ratio sheer viscosity to entropy 
density rj/s, and led to the announcement of the discovery of the nearly per- 
fect fluidity of the sQGP [1, 2]. The key role of the viscous corrections has 
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been also addressed already by the very first V2 data for Pb-Pb coUisions at 
the LHC [3]. Indeed, the shear viscosity plays a crucial role in the estimation 
of the initial condition of the heavy ion collisions [4]. 

Additionally, the hybrid model VISHNU that couples viscous hydrody- 
namics for the macroscopic expansion of the QGP to the hadron cascade 
model for the microscopic evolution of the late hadronic stage is a very 
promising tool for a more precise extraction of the QGP shear viscosity [5]. 
Therefore, several studies on shear and bulk viscosity have been performed 
in the last few years, for the confined [6, 7, 8] as well as for the deconfined 
phase [9, 10, 11, 12, 13]. Additionally, also for the thermal conductivity mod- 
els have been developed [14, 15]. Therefore a dynamic consistent calculation 
of all transport coefficients, i.e. the shear viscosity 77, the bulk viscosity ( 
and the heat conductivity k, within a model describing the strong coupling 
properties of the QGP is mandatory. Recently, we have performed an inves- 
tigation of the viscosity r] of the sQGP [9] in a dynamical way within kinetic 
theory in the ultrarelativistic limit using the virial expansion approach in- 
troduced in Ref. [16]. In this approach the choice of the interaction between 
the partons in the deconfined phase plays a crucial role to reproduce the 
thermodynamic properties of the QGP in comparison to lattice QCD. By us- 
ing a generalized classical virial expansion we have calculated the corrections 
to a single-particle partition function starting from an interaction potential 
whose parameters are fixed by thermodynamical quantities. From the same 
potential one can evaluate the transport cross section and then the sheer 
viscosity [9]. 

We have found 77/s fa 0.097 which is very close to the theoretical lower 
bound. Furthermore, for T < l.STc the ratio is in the range of the experi- 
mental estimates 0.1 — 0.3 extracted from RHIC experiments. 

In the present work we test this approach to the bulk viscosity and the 
heat conductivity. In particular, we focus on the ultrarelativistic expan- 
sion of the transport coefficients. We use a relativistic extension [17] of the 
Chapman-Enskog method [18] developed in [19] to calculate the values of 
the transport coefficients for a gas with arbitrary particle interaction. The 
specific form of the interaction is encoded in the transport cross section. 
Very recently, this method was reviewed in Ref. [12] for the shear viscos- 
ity only. The aim of the authors was to perform a quantitative comparison 
between the results of shear viscosities from the Chapman-Enskog and relax- 
ation time methods for selected test cases with specified and simplified elastic 
differential cross sections. In this contribution we use the Chapman-Enskog 
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method for one specific quark-quark interaction. This leads to an equation of 
state in fine with recent three-flavor QCD lattice data for the pressure, speed 
of sound, and interaction measure at nonzero temperatures and vanishing 
chemical potential. Furthermore, from the same interaction we have directly 
extracted the effective coupling ay to be employed for the determination of 
the transport cross section which enters the transport coefficients. In this 
way, we systematically calculate all three transport coefficients 1], ( and k,. 

This paper is organized as follows. In Section 2, following Refs. [12, 20] 
we review the formal expression for the transport coefficients in the first 
order of the first Chapman-Enskog approximation. In particular, we give a 
formulation in terms of relativistic omega integrals. In the ultrarelativistic 
limit these integrals arc evaluated for cross sections which are independent 
on the relative and total momentum of two colliding particles. 

In Section 3 we calculate and discuss the transport coefficients for our 
model within the ultrarelativistic approximation. The conclusions in Sec- 
tion 4 finalize this work. 



2. Theory 

A elegant and clear derivation of the general expression for the trans- 
port coefficient in the first Chapman-Enskog approximation can be found in 
Ref. [19]. In the following we directly recall the results for the heat conduc- 
tivity and the sheer and bulk viscosity in the first order of this approximation 
of a gas of free particles with mass m at temperature T 

3 m Oil 

C - T^. (3) 

0^22 

The quantities a2, Pi and 70 are known functions of the temperature and 
independent of the particle interaction. In general, by introducing the di- 
mensionless quantity 

z^m/T, (4) 

the enthalpy per particle 
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and the adiabatic coefficient 

7=^ (6) 
cv 



tfiese scalar function are given by 



70 = -10/1 = -10-^, (7) 

A = -3^ = d(;.3(.)/i^.(2)} 
7 — 1 d^;"^ 
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Q!2 = 7: 



^^(7-^1+7 



(9) 



wliere Kn{z) denotes tlie modified Bessel function of tlie second kind of order 
n. Tfie quantities 022, &11 and cqo are interaction independent and formally 
given by 

coo = 16 - .-^n ) , (10) 

6n = 8 ) + ) , (11) 

022 = 2a;f , (12) 
where the relativistic omega integrals are in general defined as 

uj\^\z) — Jjl^ , / du sinh^w cosh* u (22; cosh it) 
Ki{z) Jo 

X ^d^sin^^^(l-cos^^), j-l + \{-iy.{lS) 



The particle interaction is encoded in the differential cross section da{u, 6) / dfl. 
Here 6 denotes the scattering angle in the center of momentum frame. The 
variable u is directly connected to the relative momentum qi2 and total mo- 
mentum Q of the two colliding particles by 

gi2 = msinhM and Q = 2mcosh.u. (14) 

One has to evaluate the omega integrals in order to calculate the transport 
properties of the system. The results, in particular, the behavior of the cross 
section as function of u and 9, are of course model dependent. The simplest 
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case, i.e. constant cross section, which correspond to an assumption of hard 
spheres, has been investigated in Refs. [12, 20]. In following we consider 
the more general case, where the cross section is momentum independent, 
i.e. da/{dfl) = da{9)/dfl and we evaluate the omega integrals under this 
assumption. 

2.1. Omega integrals in ultrarelativistic limit 

The momentum independence of the differential cross section leads to a 
simplification in Eq.(13). Performing the the integration the u integration, 
using X — cosh rt, yields 

where the angular integral VF^*) is given by 

W^'^ = I de sine (1 - cos^^). (16) 

Jo diZ 

We note here, that for the calculation of the transport coefficients omega 
integrals for s = 2 have to be computed. Then, one can directly connect the 
integral W^'^'> with the (viscous)transport section defined by [21] 



One obtains the relation 



= fdn^ sm^e. (17) 

J dil 



27T 

and the relativistic omega integral may be written as 



= ^ (18) 

2n 



oo 



uP{z) = j^^z"^ I dx (x^ - If x' Kj{2zx). (19) 



By replacing zx ^ x and using the binomial formulae we obtain 



3 /o\ coo 



Because the limit z ^ corresponds to the ultrarelativistic limit, an ex- 

4' 



(2) 

pansion of oo) around z = leads to a systematic evaluation of the transport 



coefficients in tfiis approximation. A detailed investigation of tfiis expansion 
allows to estimate the importance of the leading order terms as well as to 
systematically compare the different transport coefficients at the same order 
in z. In fact, in the ultrarelativistic limit, the integral in Eq.(20) becomes 
independent of because the lower bound can be replaced in this expres- 
sion by zero. This replacement allows to perform the integration analytically 
using 

/ Ax K,{2x) = -r[(// + iy + l)/2]r[(/i - u + l)/2] (21) 
Jo 4 

for ± 1/ + 1 > (see Appendix A), for the calculation of the transport 
coefficients this condition is satisfied. In the ultrarelativistic limit the omega 
integral may be written as 

"^'''^^^ = iK^U) ,^.4 E (^)(-l)V^r[(;.+.+l)/2]r[(;.-.+l)/2]. (22) 

We emphasize here, that this result is a fractional expression, where the 
numerator is given by a polynomial in z. The denominator contains the 
modified Bessel function K2{z), whose expansion at 2; = can not be ex- 
pressed as a polynomial as shown in Appendix A. For K2{z) the leading 
terms of the expansion are 



Mz) = ^ 



z^ 



1 - 7^' + -4^' + ^CEM-z' + Oiz"" In z) 
4 ib 62 



(23) 



where the constant cem is connected to the gamma of Euler-Mascheroni 7em 

by 

3 

Cem = 2 In 2 + - - 27em ~ 1.731863. (24) 
We obtain for the omega integrals the the following expansion in z 

^f\^) = ^r[(z+j + 7)/2]r[(z-j + 5)/2] 

X \z^^- +o(z^ Mnz), (25) 

V i2-j2 + 10i + 25 2 y ^ ^ ' 

where only the first two leading terms are considered. Combining this ex- 
pression with the analytic expression for the quantities «2 ■, f^i and 70 given 
in Eqs.(7-9) we can analytically evaluate the transport coefficients. 
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2.2. Evaluation of transport coefficients in ultrarelativistic limit 

For the calculation of the transport coefficients in the ultrarelativistic 
limit we also need the asymptotic expansion for the scalar function «2) Pi 
and 7o. Using Eqs.(7-9), they can be expressed in term of the Bessel function 
K2{z) and -^^3(2;) In addition to the asymptotic expansion K2{z) given in 
Eq.(23), we need the corresponding expression for K^lz), 



—z^ + 0{zHnz) 
d4 



(26) 



Then we obtain for the (relevant) scalar quantities 



1600 



7o = 
2 



l + -z^]+0{zHnz), 



Pi 



144 \^--z^j +0(/ln^), 

—z^ + Oiz^hiz). 
36 ^ ' 



(27) 
(28) 
(29) 



The dynamical quantities 022,^11 and Cqo defined in Eqs.(10)-(12) can be 
evaluated in the ultrarelativistic limit using Eq.(25). We obtain 



Coo = 
&11 = 
022 = 



200at A ^ 1 2^nMi ^ 
(^1+5^ +0{z \nz^ 

^(l + 0(.^ln.)), 

^(l-i.^ + 0(.Mn.) 



(30) 
(31) 
(32) 



Finally, using Eqs.(l)-(3) we can write the transport coefficients as 



77 = 



4m 



l + ^^' + 0(/ln^)), 



1 



3(7t 

mz" 
108(7t 



{l-=-z' + 0{znnz) , 



{l + 0{zHnz)) . 



(33) 
(34) 
(35) 
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Using the definition of z for the leading term of the shear viscosity we recover 
the known expression [9, 21] 

We note that the leading term is temperature dependent, but independent 
of the mass. For the heat conductivity k the leading term does not depend 

on z. Nevertheless, a implicit temperature dependence can be encoded via 
the transport cross section. Additionally, we note that in this ultrarelativis- 
tic approach the bulk viscosity is automatically suppressed in comparison to 
the other transport coefficient as expected from general physical considera- 
tions [22]. 

3. Results 

For the final calculation using Eqs. (33-35) we have to specify an interac- 
tion. We follow Ref. [9] to calculate the from the interaction Wyi used to 
describe within a virial expansion the three-flavor thermodynamics properties 
of the QGP from the lattice [23]. We use an effective quark-quark potential 
Wyi motivated in Ref. [16] and inspired by a phenomenological model which 
includes non-perturbative effects from dimension two gluon condensates [24] . 
The effective potential between the quarks reads 

W^^{r,T) = f-- + e-^(^>, (37) 

where C2 = (0.9GeV)^ is the non-perturbative dimension two condensate and 
M(T) the Debye mass estimated as 

M(T) = ViVc/3 + 7Vf/6 gT = ~gT. (38) 

Setting the coupling parameter g = 1.30 one describes the recent three-flavor 
QCD lattice data for all thermodynamic quantities [23] in the temperature 
range up to 5 Tc very well. 

The transport cross section is given by [25, 26] 

at{s) = (To 45(1 + z) [{2z + 1) ln(l + 1/z) - 2] , (39) 

with the total cross section 0-0(5) = 9710^(3) /2fi'^,.j.. Here cty = Q;v(7") is the 
effective temperature-dependent coupling constant and and z = M{TY/s. 
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For simplicity, we will assume ctq to be energy independent and neglect its 
weak logarithmic dependence on s in the relevant energy range and set s ~ 
17T2 [27]. 

In this study we consider only elastic cross sections. An extension to 
inelastic scattering processes for an ultrarelativistic Boltzmann gas and its 
application for the calculation of the shear viscosity has been recently pre- 
sented [28]. Following Refs. [9, 29] we can extract the coupling ay from the 
interaction Wu. We define the coupling in the so-called qq-scheme, 

. _H,.1!M-I). (40) 

71 dr 




Figure 1: (Color online) The transport cross section normalized to its value at the critical 
temperature crt/cr[o) as a function of the temperature expressed in units of the critical 
temperature T^. 

The coupling aqq(r, T) then exhibits a maximum for fixed temperature at 
a certain distance denoted by rmax- By analyzing the size of the maximum 
at Tmax we fix the temperature dependent coupling ay(T) as 



ay(T) = aqq(rmax,T) . (41) 



With these results we obtain a temperature dependent transport cross section 

In Fig. 1 the transport cross section normahzed to its value at the critical 
temperature, is shown as a function of the temperature expressed in units of 
the critical temperature Tc. We note a strong decrease of the cross section 
with temperature. This can be not neglected in a consistent evaluation of 
the transport coefficients. 




T/T 



Figure 2: (Color online) The scaled sheer viscosity rj/T^ (red line), heat conductivity k/T^ 
(blue line) and bulk viscosity C/T^ (green line) as function of the temperature expressed 
in units of the critical temperature Tc . 

In order to investigate the importance of the different coefficients we 
scale them by the appropriate power of the temperature in order to obtain 
a dimensionless quantity. Therefore we show in Fig. 2 the ratios rj/T^ (red 
line), k/T^ (blue line) and (/T^ (green line) as a function of the temperature 
expressed in units of the critical temperature T^. Several features become 
evident: 

i) For the shear viscosity and the heat conductivity an increase with the 
temperature is found for Tc < T < 3.5Tc. 
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Figure 3: (Color online) The bulk viscosity ( (solid line) as function of the temperature 
expressed in units of the critical temperature Tc. In comparison the bulk viscosity using 
the temperature independent transport cross section ctj"^ . 

ii) At higher temperatures ?7/T'^ and k/T^ become constant. This behavior 
of the shear viscosity is also present in very different approaches, such 
as AdS/CFT [30], the quasiparticle approximation including the quark 
selfenergy [31, 32] and in the weak coupling estimate from Ref. [33, 34]. 
For the heat conductivity that is also found in perturbative calcula- 
tions [14, 35], where k, scales with the temperature square. 

Hi) The bulk viscosity is completely negligible in comparison to the other 
coefficients. This justifies the omission of ( in several hydrodynamical 
and transport calculations [36, 37, 38]. 

In Fig. 3 we show ( as function of the temperature expressed in units of 
the critical temperature Tc (solid green line). To investigate the role of the 
temperature dependent transport cross section we compare with the results 
obtained using the constant cross section (t|°'' (green dashed line) In the 
calculation with the full cross section the suspected enhancement of the bulk 
viscosity is missing. In contrast, the interplay between cross section and 
leads to a broad maximum around 3.2Tc. This is clearly a consequence of 
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the temperature dependence of Ut: namely, using the constant value a^'^^ for 
the transport cross section, we find a decreasing behavior, that suggests the 
possibility of a maximum of the bulk viscosity at T^. It is not surprising 
that our ultrarelativistic model can not describe the peak of ( at the critical 
temperature. This maximum is a consequence of the hadronic correlations at 
Tc [39] that are non included in this approach. To include such correction one 
can add to this parameter free derivation of the cross section an estimation 
of the hadronic transport cross section and the omega integrals for the whole 
temperature range have to be calculated. Therefore, analytic expressions for 
the coefficients can not performed [20]. 

Alternatively, an estimate for the correlation term of the the bulk viscosity 
can be given using the expression (cov = ~A{Cs — 1/3)^ , with A = 4.558 or 
A = 4.935 [40] or A = 2 [41] respectively. 

Another important aspect of our results regarding the sheer viscosity and 
the heat conductivity is the relative importance of the next to leading order 
contributions. In Fig. 4 the percent deviation Arj of the sheer viscosity (red 



0.02 



^ 0.015 




T/T 



Figure 4: (Color online) The percentage deviation sheer viscosity Arj (red line) and of the 
heat conductivity An (blue line) as function of the temperature expressed in units of the 
critical temperature Tc. 

line) and A/t of the heat conductivity (blue line) have been investigated. The 
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percent deviations are defined by 



AX = J ^ X 100, (42) 

X 

where X = rj, k and denotes the next to leading order contribution. 

Evidently, the deviations from the leading order terms are very small and 
the higher order contributions can be neglected. In this sense, the calcula- 
tions of the ratio sheer viscosity to entropy density performed in Ref. [9], 
where the leading term of the uhrarclativistic expansion and the entropy 
within this virial expansion model have been used, is unchanged by the higher 
term corrections. 



4. Conclusions 

We presented the calculation of the transport coefficients of the QGP by 
revisiting the first order of the first Chapman- Enskog approximation. We 
formulated them in terms of the relativistic omega integrals. Their evalua- 
tion in the ultrarelativistic hmit is exphctly shown by assuming that the cross 
section is independent on the relative and total momentum of two colliding 
particles. The well known features of the ultrarelativistic hmit are recov- 
ered, i.e. the suppression of the bulk viscosity in comparison to the other 
coefficients. In particular, the expected peak of C at the critical temperature 
is missing, because the hadronic correlations at have not been included. 
However, the overall suppression justifies neglecting of the bulk viscosity in 
several calculations for the evolution of the QGP [4, 36, 37, 38]. For the sheer 
viscosity and heat conductivity we observe an increasing behavior with the 
temperature. At high temperatures, they scale as power law in T. Addition- 
ally, we find that all results in this ultrarelativistic expansion are dominated 
by the leading order contribution. Therfore, for the ratio sheer viscosity to 
entropy density we recover the results of our previous estimation [9] . 

For further work, the limitation of the ultrarelativistic limit has to be re- 
laxed and the coefficients can be evaluated in the whole temperature range. 
This allows the including of hadronic correlations at in the calculations of 
(. Additionally, we can implement the temperature dependent transport co- 
efficients calculated here as well as our realistic equation of state in dissipative 
hydrodynamical calculation. In this way we can systematically investigate 
the evolution of the QGP in the heavy-ion collisions following the lines of 
Ref. [4]. 
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A. Modified Bessel functions of second kind 

The Bessel functions [42], first defined by the mathematician Daniel 
Bernoulli and generalized by Friedrich Bessel, are canonical solutions w{z) 
of Bessel's differential equation 

^2 -5-^ + X— + ^ - w = 43 
dz^ dz 

for an arbitrary real or complex number v which indicate the order of the 
Bessel function. In general, v is integer or half-integer. The solutions called 

Bessel functions of the first kind are denoted as Jyi^z) and arc holomorphic 
functions of z throughout the 2;-planc cut along the negative real axis. It is 
possible to define the function by its Taylor series expansion around 2; = 

where V{z) is the Gamma function defined by 

T{z) = / t'-^e-'dt. (45) 



By replacing z by ±iz one obtains the modified Bessel equation 

z^-r^ + X— + {z^ + iy^)w = 0. 46) 
dz'^ dz 

The solutions are called modified Bessel functions of first and second kind, 
and K^{z) respectively, and are connected to the functions Ju{z) by 

Uz) = t-'-Mtz), (47) 
= |M£l_^ (48) 

Unlike the ordinary Bessel functions, which are oscillating as functions of a 
real argument, the modified Bessel functions are exponentially growing and 
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decaying functions, respectively. For real arguments the function Iu{x) 
goes to zero at a; = for z/ > and is finite at a; = for z/ = 0, like the 
ordinary Bessel function Ju{x). Analogously, Ki,{x) diverges at x = 0. 

There are many integral representations of these functions, using integrals 
along the real axis as well as contour integrals in the complex plane. In the 
following we focus on the modified Bessel functions of the second kind. Very 
useful integral representations are 

/ZZ/l^\u poo -| 1 

Ku{z) = Y^^^Jdte-\t'-ir-K 5^^>-2, |ph ;.| < -tt (49) 

POO -r 

K^{z) = / dt e-^"°^^*cosh(i/t), |ph;2|<-7r. (50) 
Jo 2 

Like the ordinary Bessel function, the modified Bessel function of second kind 
can be expressed as a power series. In the case of integer order the expansion 
reads 



fc=0 



2'2 ' ^ A;! ^ 4 ^ ^ ^ ^2 



where ip{w) = T'{w)/T{w) denotes the ip function, which satisfies 

n ^ 

t/^(l) = -7em and ^'(n + 1) = ^ - - 7em- (52) 

fc=i 

Here 7em denotes the Euler-Mascheroni constant. The expansion for K2{z) 
and Ks{z) is obtained by combination of Eq. (44) and (47) together with 
the definitions in Eq. (52). In addition to the normalization-like condition 



I 



oo 11 

dt K^{t) = -7rsec(-7ri/), \'Rv\<l (53) 



we recall here 



which is equivalent to the expression in Eq. (21). 
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